NBI-HE-98-11 
AEI-064 
MPS-RR-98-5 



00 
On 



OO 



Oh 



X 



Non-perturbative Lorentzian Quantum Gravity, 
Causality and Topology Change 



J. Amhj0rn a >f\ and R. Lo]J 6 f] 



& • a The Niels Bohr Institute, 



Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark, 

6 Max-Planck-Institut fur Gravitationsphysik, 
Albert-Einstein-Institut, 



00 ' Schlaatzweg 1, D- 14473 Potsdam, Germany 

O ' 

in 
o 

00 
On 



Abstract 



We formulate a non-perturbative lattice model of two-dimensional Lorentzian quan- 
tum gravity by performing the path integral over geometries with a causal structure. 
The model can be solved exactly at the discretized level. Its continuum limit coin- 
cides with the theory obtained by quantizing 2d continuum gravity in proper-time 
gauge, but it disagrees with 2d gravity defined via matrix models or Liouville the- 
ory. By allowing topology change of the compact spatial slices (i.e. baby universe 
creation), one obtains agreement with the matrix models and Liouville theory. 
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1 Introduction 



A moot point in non-perturbative Euclidean path-integral approaches to quantum 
gravity is the final interpretation of their results in terms of physical quantities 
defined in the Lorentzian sector of the theory. Even classically, it is known that 
the simple expedient of applying a Wick rotation t — > r = it fails in all but a few 
special cases, for instance, when the space-time is static and thus admits a global 
choice of time Xq such that the cross terms goi, i > 1 of the metric tensor vanish 
and its spatial components g^ are time-independent. We do not know of a way to 
set up a 1-to-l correspondence between generic solutions of the continuum Einstein 
equations with different signatures. If it exists, it is likely to be technically involved 
(see, for example, QIJ for a recent proposal). 

From this point of view it is unclear if one can expect to obtain the correct 
Lorentzian theory by first performing a path integral over general Euclidean metric 
configurations and then analytically continuing in some way. Related questions have 
been discussed in 0, where it was suggested that in the path integral for gravity 
with an action in square-root form, and using the rather unconventional weights 
e^ 5 , it may be necessary to sum over both Lorentzian and Euclidean metrics in 
order to obtain a unitary evolution. 

To investigate these issues further, we formulate a theory of gravity in two space- 
time dimensions where the summation is restricted to metric configurations with a 
causal structure, defined on each contribution to the sum over states in the discrete 
model under consideration. In this way one encodes at least part of the Lorentzian 
structure into the Euclidean path integral. In addition to that, we will consider a 
suitable analytic continuation of our results.^ Curiously, our final continuum results 
show some similarity with the approach advocated in []2| , although our path-integral 
construction proceeds along entirely different lines. 

The idea that a notion of causality should be built into each history that con- 
tributes to the path integral amplitude goes back at least to Teitelboim [[|, and 
has more recently been advocated in || ||. However, to our knowledge a concrete 
implementation in a well-defined, non-perturbative model for quantum gravity has 
so far been missing. An ideal testing ground for this idea is in two dimensions, 
where the Euclidean path-integral construction leads to a non-trivial gravitational 
quantum theory. Its properties have been explored in great detail, both by ana- 
lytical and numerical methods. The main reason why this model can be solved 
analytically, even at the discretized level, is that the action in 2d gravity is trivial. 
The Einstein-Hilbert term is a topological invariant and does not contribute unless 
we consider topology changes of space-time. Moreover, for fixed space-time volume 
the partition function is purely entropic and given by the number of different geo- 
metries. Exactly this fact is used in the formalism of dynamical triangulations (or 

3 Another way of making 2d gravity more Lorentzian (by allowing for both time- and spacelike 
edges in a dynamically triangulated model, but without causal structures) was studied in ||, but 
no significantly different behaviour from the Euclidean theory was found. 
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equivalently, matrix models) to construct the non-perturbative path integral. One 
counts the number of inequivalent triangulations which can be constructed from a 
given number of triangles of unit volume, and with a fixed topology for the result- 
ing simplicial complex. Inequivalent triangulations (appropriately defined) can be 
related to different geometries |7|, H [HJ. In this way the following quantities have 
been calculated in pure Euclidean 2d gravity: 

(1) the partition function on the sphere, (A), as a function of the renormalized 
cosmological|7| constant A, 

Z {eu) (K) ~ A 2-7 + terms less singular in A, 7 = --, (1) 

2 

where 7 is the so-called string susceptibility exponent; 

(2) the Hartle-Hawking wave functional W\(L) as a function of the length L of 
the spatial boundary (with a marked point) Hill , 



Wt\L)Kj^(l + VAL)e-^ L , (2) 
(more generally, one can calculate "multi-loop" correlators W^(Li, . . . ,L n ) 

0); 

(3) the average Euclidean space-time volume, By(R), of geodesic balls of radius 



R in the ensemble of universes of Euclidean space-time volume V [13, 14, 15 



B v (R)kR*f(^), (3) 

where F(0) = 1 and F(x) oc e~ xl/3 for x — > 00. Relation @ implies that the 
intrinsic Hausdorff dimension dn of Euclidean 2d quantum gravity is four (and 
not two, as one naively might have expected) .0 

In the following we will show that for a universe with cylinder topology, restrict- 
ing the path integral to configurations admitting a causal structure leads to a theory 
with 7 = 1/2 (although the definition of 7 turns out to be ambiguous), a Hartle- 
Hawking wave function with the same exponential decay (in the Euclidean sector) 
as (0), but a different functional form, and an intrinsic Hausdorff dimension of two, 
and not four as in (|3|) . We will further show that once we allow for topology changes 
of space, i.e. the creation of baby universes, we are led to a theory satisfying (l)-(3). 

4 It can be shown, using only general arguments ]T^ | , that Euclidean space-time, with an assumed 
fractal dimension du, is characterized by a function like (0): 

B V (R) - R dH F[R dH /y), F(0) = 1, F(x) = e^"^ 1 ' for x > 1. (4) 
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Figure 1: The propagation of a spatial slice from step t to step t + 1. The ends of 
the strip should be joined to form a band with topology S 1 x [0, 1]. 

2 The discrete model 

As mentioned in the introduction, the solution of two-dimensional quantum gravity 
amounts to counting geometries. While this counting problem has been solved in 
Euclidean gravity, it seems non-trivial if the space-time has Lorentzian signature. 
Counting in a field theoretical context usually amounts to the introduction of a 
regularization (a discretization) which makes the counting procedure well defined. 
After the counting has been performed one may attempt to take the continuum limit 
of the discretized theory. It is unclear which class of geometries to include if the 
signature is not Euclidean. We propose here a model where a causal structure is 
explicitly present in all the geometries included in the path integral. 

The model is defined as follows. The topology of the underlying manifold is taken 
to be S 1 x [0, 1], with "space" represented by the closed manifold S 1 . We consider 
the evolution of this space in "time" . No topology change of space is allowed at this 
stage, but we will return to this issue in sec. |5|. 

The geometry of each spatial slice is uniquely characterized by the length assigned 
to it. In the discretized version, the length L will be quantized in units of a lattice 
spacing a, i.e. L = I ■ a where I is an integer. A slice will thus be defined by I 
vertices and I links connecting them. To obtain a 2d geometry, we will evolve this 
spatial loop in discrete steps. This leads to a preferred notion of (discrete) "time" 
t, where each loop represents a slice of constant t. The propagation from time-slice 
t to time-slice t + 1 is governed by the following rule: each vertex % at time t is 
connected to fcj vertices at time t + 1, hi > 1, by links which are assigned length —a. 
The hi vertices, hi > 1, at time-slice t + 1 will be connected by fc« — 1 consecutive 
space-like links, thus forming hi — 1 triangles. Finally the right boundary vertex in 
the set of hi vertices will be identified with the left boundary vertex of the set of 
k i+1 vertices. In this way we get a total of Y?i=i{h — 1) vertices (and also links) 
at time-slice t + 1 and the two spatial slices are connected by J2\=i h = k + k+i 
triangles. See fig. [1]. 

The elementary building blocks of a geometry are therefore triangles with one 
space- and two time-like edges. We define them to be flat in the interior. A consistent 



way of assigning interior angles to such Minkowskian triangles is described in [|17 
The angle between two time-like edges is j t t — ~ arccosf, and between a space- 
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and a time-like edge j st = f + \ arccos |, summing up to j t t + ^Ist = 7T- The sum 
over all angles around a vertex with j incoming and k outgoing time-like edges (by 
definition j, k > 1) is given by 2n + (4 — j — k) arccos |. The regular triangulation 
of flat Minkowski space corresponds to j = k = 2 at all vertices. The volume of a 
single triangle is given by ^a 2 . 

One may view these geometries as a subclass of all possible triangulations that 
allow for the introduction of a causal structure. Namely, if we think of all time-like 
links as being future-directed, a vertex v' lies in the future of a vertex v iff there is 
an oriented sequence of time-like links leading from v to v'. Two arbitrary vertices 
may or may not be causally related in this way. 

In quantum gravity we are instructed to sum over all geometries connecting, say, 
two spatial boundaries of length L\ and L 2 , with the weight of each geometry g 
given by 

e^M S[g) = AjV=9 (in 2d), (5) 

where A is the cosmological constant. In our discretized model the boundaries 
will be characterized by integers l\ and l 2 , the number of vertices or links at the 
two boundaries. The path integral amplitude for the propagation from geometry 
li to l 2 will be the sum over all interpolating surfaces of the kind described above, 
with a weight given by the discretized version of @. Let us call the corresponding 
amplitude G^\li,l 2 )- Thus we have 

00 

G\\h, h) = J2 G x\h,k;t), (6) 



00 



G^(h,l 2 ;t) = £Gi 1) (Ji,/;l)/G? ) (J,M-l), (7) 
1=1 



G 



(h,l 2 ;l) = r E e^E^, (8) 



(i) 

x \ ■ 1 ■ ■ . • ■ 1 j 

1 {ki,...,k h } 



where A denotes the bare cosmological constant^ (we have absorbed the finite triangle 
volume factor), and where t denotes the total number of time-slices connecting l\ 
and l 2 - 

From a combinatorial point of view it is convenient to mark a vertex on the 
entrance loop in order to get rid of the factors / and l/l in (|7|) and (|8|), that is, 

G x {hMt)=kGf{h^t) (9) 

(the unmarking of a point may be thought of as the factoring out by (discrete) 
spatial diffeomorphisms) . Note that G\{h, h]^) plays the role of a transfer matrix, 

5 One obtains the renormalized (continuum) cosmological constant A in (^|) by an additive 
renormalization, see below. 
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satisfying 

Gxiluk^ + h) = E G A(/i,/;ti)G A (/,/ 2 ;t 2 ) (10) 

i 

G x (h,k;t + 1) = J2 G *(h,l;l)G x (l,l 2 ;t). (11) 

i 

Knowing Gx(h,h'A) allows us to find G^l^^t) by iterating ([TT|) t times. This 
program is conveniently carried out by introducing the generating function for the 
numbers Ga(^i, h] t), 

G x (x,y;t)=J2x k y l G x (k,l;t), (12) 

k,l 

which we can use to rewrite fllCf ) as 

/dz 
^— G x (x, z~ ; ti)G x (z, y; t 2 ), (13) 

where the contour should be chosen to include the singularities in the complex bi- 
plane of G x (x, z~ x \ ti) but not those of G x (z,y;t2). 

One can either view the introduction of G x (x,y;t) as a purely technical device 
or take x and y as boundary cosmological constants, 

x = e iX * a , y = e lX ° a , (14) 

such that x k = e lXiak becomes a boundary cosmological term, and similarly for 
yi = e i\ ai_ us f or notational convenience define 

g = e iXa \ (15) 

For the technical purpose of counting we view x, y and g as variables in the complex 
plane. In general the function 

G(x,y;g;t) = G x (x,y;t) (16) 

will be analytic in a neighbourhood of (x, y, g) = (0, 0, 0). 

From the definitions (^|) and @ it follows by standard techniques of generating 
functions that we may associate a factor g with each triangle, a factor x with each 
vertex on the entrance loop and a factor y with each vertex on the exit loop, leading 
to 

oo / oo \ ^ oo 2 

°<™« 1] ^ S [«*£<">') - £<«*>' = (17) 

Formula (0) is simply a book-keeping device for all possible ways of evolving from 
an entrance loop of any length in one step to an exit loop of any length. The 
subtraction of the term 1/(1 — gx) has been performed to exclude the degenerate 
cases where either the entrance or the exit loop is of length zero. 
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From ([LTD and eq. (|13"D, with t\ = 1, we obtain 

G(x, y; g; t) = G(-°—, y;g;t-l). (18) 

1 — gx 1 — gx 

This equation can be iterated and the solution written as 

2 

G(x,y;g;t) = F*(x)F*(x) ■ ■■F t 2 _ 1 (x)- f %V (19) 

where F t (x) is defined iteratively by 

F <W = i- g FUxY F " (x) = x - m 

Let F denote the fixed point of this iterative equation. By standard techniques one 
readily obtains 

t[X) ~ 1- X F + F^(x-Fy 2g " 1 ] 

Inserting (^T|) in eq. we can write 

F 2 \\-F 2 ) 2 xy 

(X ' m 9 " t] ~ (A t -B t x)(A t -B t (x + y) + C t xy) (22) 



F 2t (l - F 2 ) 2 xy 



(23) 



(l-xF)-F 2t + 1 (F-x)\ [(l-xF)(l-yF)-F 2t (F-x)(F-y) 

where the time-dependent coefficients are given by 

A t = l-F 2t+2 , B t = F(l-F 2t ), C t = F 2 (l- F 2t ~ 2 ). (24) 

The combined region of convergence to the expansion in powers g k x l y m , valid for all 
t is 

\g\ < W < \v\ < L ( 25 ) 

The asymmetry between x and y in the expressions ( |2"2"|) and ( |2"3"D is due to the 
marking of the entrance loop. If we also mark the exit loop we have to multiply 
G\(h,h',t) by l 2 - We define 

G%\h, l 2 - 1) = l 2 G x (h, l 2 - 1) = hl 2 G^{h, l 2 - 1). (26) 

The corresponding generating function G^(x,y; g;t) is obtained from G(x,y; g;t) 



by acting with y-p, 



g*>, x „. e . t) = f (i-f)'*y (27) 



We can compute Gx{h,h',t) from G(x,y; g;t) by a (discrete) inverse Laplace 
transformation 



dx r dy 11 
2ni x J 2%i y x l1 y l 



G *(h, k; t)=f ^— fTdh.-^ — G & v\ 9; t), (2? 



where the contours should be chosen in the region where G(x, y; g; t) is analytic. A 
more straightforward method is to rewrite the right-hand side of ( p2|) as a power 
series in x and y, yielding 



k 



F 2 *(1-F 2 ) 2 E^ ntaMH i 1+ i 2 - k -l ( A t C t . 

Al l ' 2, j i 2 M(h-k-m 2 -k-i)\ [ b? ) ■ 

(29) 

which, as expected, is symmetric with respect to l\ and I2 after division by l\. 

In the next section we will give explicit expressions for G\(h, h',t), G\{li, I2) and 
G\(x,y) (the integral of G\(x,y;t) over t) in a certain continuum limit. 



3 The continuum limit 

The path integral formalism we are using here is very similar to the one used to re- 
present the free particle as a sum over paths. Also there one performs a summation 
over geometric objects (the paths), and the path integral itself serves as the propa- 
gator. From the particle case it is known that the bare mass undergoes an additive 
renormalization (even for the free particle), and that the bare propagator is subject 
to a wave-function renormalization (see f\M for a review). The same is true in two- 



dimensional gravity, treated in the formalism of dynamical triangulations The 
coupling constants with positive mass dimension, i.e. the cosmological constant and 
the boundary cosmological constants, undergo an additive renormalization, while 
the partition function itself (i.e. the Hartle-Hawking-like wave functions) undergoes 
a multiplicative wave-function renormalization. We therefore expect the bare cou- 
pling constants A, Aj and Ao to behave as 

A = % + A, A, = ^ + X, \ = ^l + Y, (30) 
a z a a 

where A, X, Y denote the renormalized cosmological and boundary cosmological 
constants. If we introduce the notation 

g c = e lC \ x c = e iC \ y c = e iC >», (31) 

for critical values of the coupling constants, it follows from (0) and ([15]) that 



g = g c e ia2K , x = x c e iaX , y = y c e iaY . (32) 
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Figure 2: The circle of convergence in the complex g plane (radius 1/2), and the 
critical lines, ending in g = ±1/2. 



The wave-function renormalization will appear as a multiplicative cut-off dependent 
factor in front of the bare "Green's function" G(x, y\ g; t), 

G A (X,Y;T) = \imy G(x,y;g;t), (33) 

where T = at, and where the critical exponent rj should be chosen so that the 
right-hand side of eq. (^) exists. In general this will only be possible for particular 
choices of g c , x c and y c in fl3"Bp. 

The basic relation (|T0|) can survive the limit (|33|) only if rj — 1, since we have 
assumed that the boundary lengths L\ and L 2 have canonical dimensions and satisfy 

Li = a li . 

From eqs. (^) and (^4|) it is clear that we can only obtain a non-trivial continuum 
limit if \F\ — > 1. This leads to a one-parameter family of possible choices 

g c = — - — for F = e ia , a G H, (34) 
2 cos a 



for critical values of g. It follows from (1TB) that most values of g c correspond to 



a complex bare cosmological constant A. However, the renormalized cosmological 
constant A in ( PU[ ) (depending on how we approach g c in the complex plane) could 
in principle still be real. 

A closer analysis reveals that only at g c = ±1/2, corresponding to a = 0,7r, is 
there any possibility of obtaining an interesting continuum limit. Note that these 
two values are the only ones which can be reached from a region of convergence of 
G(x,y; g;t) (see fig. 0). Note also that requiring the bare A to lie inside the region 
of convergence when g — > g c leads to a restriction ImA > on the renormalized 
cosmological constant A, since \g\ < | =^ Im A > In 2. 

Without loss of generality, we will consider the critical value g c = 1/2. It corre- 
sponds to a purely imaginary bare cosmological constant A c := C\/a 2 = —i\n2/a 2 . 
If we want to approach this point from the region in the complex g-plane where 
G(x, y; g; t) converges it is natural to choose the renormalized coupling A imaginary 



as well, A = zA, i.e. 



ln± 

A = % ± iA. (35) 



a 2 
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One obtains a well-defined scaling limit (corresponding to A 6 1R) by letting A — > A c 
along the imaginary axis. The Lorentzian form for the continuum propagator is 
obtained by an analytic continuation A — > — iA in the renormalized coupling of the 
resulting Euclidean expressions. 

At this stage it may seem that we are surreptitiously reverting to a fully Euclidean 
model. We could of course equivalently have conducted the entire discussion up to 
this point in the "Euclidean sector" , by omitting the factor of —i in the exponential 
(|) of the action, choosing A positive real and taking all edge lengths equal to 1. 
However, from a purely Euclidean point of view there would not have been any 
reason for restricting the state sum to a subclass of geometries admitting a causal 
structure. The associated preferred notion of a discrete "time" allows us to define 
an "analytic continuation in time" . Because of the simple form of the action in two 
dimensions, the rotation 

J dx dt^-g lor -> z J dx dt euy /g^ (36) 

to Euclidean metrics in our model is equivalent to the analytic continuation of the 
cosmological constant A. 

From ( p2|) or ( p3[) it follows that we can only get macroscopic loops in the limit 
a — > if we simultaneously take x, y — > 1. (For g c = —1/2, one needs to take x, y — > 
— 1. The continuum expressions one obtains are identical to those for g c = 1/2.) 
Again the critical points correspond to purely imaginary bare boundary cosmological 
coupling constants. We will allow for such imaginary couplings and thus approach 
the critical point Aj = A D = from the region of convergence of G(x, y; g; t), i.e. via 
real, positive X, Y where 

Aj = iXa, X = iYa. (37) 

Again X and Y have an obvious interpretation as positive boundary cosmological 
constants in a Euclidean theory, which may be analytically continued to imaginary 
values to reach the Lorentzian sector. 
Summarizing, we have 

9 = \e~ Aa2 - ~(1 - ha 2 ), (i.e. F = 1 - aVX) (38) 

as well as 

x = e~ Xa -^1-aX, y = e~ aY -> 1 - aY, (39) 

where the arrows — > in (^) and ( |3~9"1) should be viewed as analytic coupling constant 
redefinitions of A,X and Y, which we have performed to get rid of factors of 1/2 
etc. in the formulas below. With the definitions (|38|) and (|39D it is straightforward 
to perform the continuum limit of G(x, y; g, t) as (x, y, g) — > (x c , y c , g c ) = (1, 1, 1/2), 
yielding 



4A e ~ 2VXT 

G A (X,Y;T) = —= — -= = 
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X (v / A + X)(v / A + V) -e- 2 ^ T {VA-X){y/X-Y)' ^ 

For T —>■ oo one finds 

A A p -2v / AT 

G A (X,F;T) T -=3? ^ (41) 

(x + Va) 2 (f + Va) 

From Ga(-^, F; 7") we can finally calculate Ga(Li, L2; T), the continuum ampli- 
tude for propagation from a loop of length L\, with one marked point, at time-slice 
T = to a loop of length L 2 at time-slice T, by an inverse Laplace transformation, 

/ioo rioc 
dX dYe XLl e YL2 G A (X,Y;T). (42) 
-ioo J—ioo 

This transformation can be viewed as the limit of (]28|) for a — » 0. The continuum 
version of (0) thus reads 

/ioo 
rfZ G A (X, -Z; TO G A (Z, F; T 2 ), (43) 
-ioo 

where it is understood that the complex contour of integration should be chosen to 
the left of singularities of G\(X, —Z\T\), but to the right of those of Ga(Z,Y,T 2 ). 
By an inverse Laplace transformation we get in the limit T — > 00 

G A (L l7 L 2 ; T) ? 4L ie -^ Ll+L2 ) e" 2 ^, (44) 

where the origin of the factor L\ is the marking of a point in the entrance loop. For 
T — > we obtain 

G A (X,F;T)^^1 7 , (45) 

in agreement with the expectation that the inverse Laplace transform should behave 
like 

G K {L 1: L 2 ;T)^ 5{L X - L 2 ). (46) 

The general expression for Ga(L\, L 2 ; T) can be computed as the inverse Laplace 
transform of formula (f40D, yielding 

GA(Ll ' L2;T) = sIuWXt" U~ h UwXtJ ' (47) 

where l\ (x) is a modified Bessel function of the first kind. The asymmetry between 
L\ and L 2 arises because the entrance loop has a marked point, whereas the exit 
loop has not. The amplitude with both loops marked is obtained by multiplying 
with L 2 , while the amplitude with no marked loops is obtained after dividing ( f4T|) 
by L\. Quite remarkably, our highly non-trivial expression ( |4~7|) agrees with the loop 
propagator obtained from a bona-fide continuum calculation in proper-time gauge 
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of pure 2d gravity by Nakayama |Tj|. More precisely, his propagator A m for m = 



(m is a winding number introduced in a somewhat ad-hoc manner in [|Tg]) is related 
to ours by Ga(Li, L 2 ; T) = -j±A (Li, L 2 ; T), which just reflects the fact that the exit 
instead of the entrance loop has been marked. Also additional ambiguities in the 
continuum formulation, involving shifts m — > m + 1/2 due to renormalization, are 
fixed in our approach. 

To obtain the propagator of the Lorentzian theory, we substitute A — > —iA 
in (|47D. As a consequence, the amplitude becomes complex and the hyperbolic 
functions pick up oscillatory contributions. Both the real and the imaginary parts 
continue to be exponentially damped for large T. What is at first puzzling about 
the functional form of (|47|) is that the naive analytical continuation in "time" , T — > 
—iT, leads to a drastically different (and highly singular) result. However, this is 
an incorrect choice, which can be understood as follows. The combination y/AT 
appearing as arguments in (O) arises in taking the continuum limit of powers of the 
form F l in expressions like fl22"D , (p3|) , where F is defined in ([H]). 

There are two aspects to a possible analytic continuation of F l . The power t in 
F l should clearly not be continued, since it is simply an integer counting the number 
of iterations of the transfer matrix. However, the function F itself does refer to the 
action, because the dimensionless coupling constant g = e tXatai is the action for a 
single Lorentzian triangle. (For added clarity we have distinguished between the 
lattice spacings in time- and space-directions, and called them a t and c^.) From 
the expression for F in terms of g in (|21~D , we have F — 1 — \fa t aiA. The analytic 
continuation of F in time, from Euclidean to Lorentzian time, corresponds to the 
substitution a t — > —i a t under the square-root sign, and thus becomes equivalent 
to the continuation A — > — iA in the cosmological constant, as already remarked 
below eq. fl36|) . The subtleties associated with the analytical continuation in the 
"time" -parameter T appearing in a transfer-matrix formulation of quantum gravity 
were first discussed in 0, |19| in the context of a square-root action formulation. 
They will be present also in more complicated theories, where the analytic continu- 
ation from Euclidean metrics to Lorentzian metrics cannot be absorbed by a similar 
continuation in A ||. 

Finally, we compute the amplitude describing the transition from L\ to L 2 for an 
arbitrary "time" -separation of the slices by integrating over T. From (f47|) , multiplied 
by L 2 in order to arrive at the symmetric propagator where both loops are marked, 
we get 

_ roc , n . p -VA|Xi-i 2 | _ p-^(L!+X 2 ) 

G ( i\L 1) L 2 )= dTGf\L u L 2 -T) = -= . (48) 

Jo 2vA 

One of course obtains the same result by first integrating G\(X, Y\ T) with respect 
to T, and then doing the inverse Laplace transform. Again the analytic continuation 
A — > —iA leads to a complex amplitude. 
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4 The differential equation 

The basic result ( pLOl) for G\(X, Y; T) can be derived by taking the continuum limit 
of the recursion relation (0). By inserting fl38|) and (^) in eq. ([18]) and expanding 
to first order in the lattice spacing a we obtain 

^G A (X, Y;T) + ^L [(x 2 - A)G A (X, F; T)J = 0. (49) 

This is a standard first order partial differential equation which should be solved 
with the boundary condition (flop at T = 0, since this expresses the natural condition 
(fl6|) on Ga(Li, L 2 ). The solution is thus 

C, (X Y- T) - ^ X )~ A 1 m 

Ga(X, Y, T) - x2 _ A x(T;X) + Y' (5 ° } 

where X(T; X) is the solution to the characteristic equation 

^- = -(X 2 -A), X(T = 0)=X. (51) 
It is readily seen that the solution is indeed given by ([!(]) since we obtain 

X(D = VA<f + *)-e-^-*), (52) 
V 7 (v / A + X) + e- 2 v / Ar( v ^ 1 _x) V ; 

If we interpret the propagator G\(Li, L 2 ; T) as the matrix element between two 
boundary states of a Hamiltonian evolution in "time" T, 

G A (L 1 ,L 2] T)=<L 1 \e- AT \L 2 > (53) 

we can, after an inverse Laplace transformation, read off the functional form of the 
Hamiltonian operator H from (j49|) . 



Using ([48]), it is now straightforward to check that 



H(L U 4 2) (L 1; L 2 ) = 5(L X - L 2 ), (55) 

as is expected for the propagator. The corresponding Hamiltonian for the propagator 
of unmarked loops is given by 

^< l 4 )= - l sp-4 +al ' < 56 > 
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A solution to the "Wheeler-DeWitt equation" is only obtained if we integrate 
the expression (|47|) over the entire T-axis (as observed a long time ago in [|J): 

oLi J-oo oL\ vA 

The above construction refers to the evolution of the system with respect to the 
"time" -parameter T appearing in the transfer-matrix approach. However, we have 
argued earlier that one should not simply analytically continue T — > — iT to relate 
the Euclidean and Lorentzian sectors of the theory. Not taking T seriously as a 
time-parameter presumably implies that also the operator H appearing in (|53|) is 
not the physically relevant Hamiltonian. 

At any rate, our choice of analytic continuation does not seem to lead to a self- 
adjoint Hamiltonian if one uses the prescription fl5"3|) for the Lorentzian case. A 
possible way out may be to use weights of the form e^ 5 and to sum over a class 
of both Lorentzian and Euclidean geometries, as advocated in [0. Like in this 
approach, we are also encountering the factor v% but it is presently unclear to us 
whether the two can be related. 



5 Observables 

We will now compare the predictions coming from the Euclidean sector of our model 
with those from 2d quantum gravity as defined by matrix models (or Liouville the- 
ory). As mentioned in the introduction, this amounts to calculating the string sus- 
ceptibility exponent 7, the Hartle-Hawking wave function W\(L) and the intrinsic 
Hausdorff dimension du- We must first define what we mean by the disc amplitude 
W\(L) in our model. A natural definition is given by 

W A (L) : =G A (L,L 2 = 0) =e-^ L , (58) 

where the last equality follows from fl48]) . We have contracted the exit loop to length 
zero in order to produce a disc. This procedure leaves a "mark" at the very endpoint 
of the universe, contrary to the usual definition of W A (L) in Euclidean 2d quantum 
gravity. Had we applied a definition analogous to (|58D in Euclidean gravity, setting 
L 2 = in the propagator would also have led to a marking inside the disc. To relate 
it to the usual disc amplitude PV{ e ^(L), formula ([58]) would have had to be replaced 
by 

- ^wt\ L ) = G { r\L,L 2 = 0) (oc T^e"^). (59) 

The derivative appears because marking a point in Euclidean gravity is equivalent 
to differentiating with respect to the cosmo logical constant. However, in our model 
the mark in (|58"D does not correspond to a differentiation d/dA, since unlike for 
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Euclidean gravity the mark cannot be located anywhere in the bulk. In fact, from 
the definition of G\(Li, L 2 ; T) it follows that 

/•OO 

W A (L)= dTG A (L,L 2 = 0;T), (60) 
Jo 

which implies that the mark is always at the latest proper time T. Although quite 
similar, (p8|) differs from ([59]), a fact which becomes more obvious when we consider 
their respective Laplace transforms, whose regular parts are given by 

w A (x) = y^-Jx> wt\x) = (X- lVK)y/x + VK. (61) 

From w[ eu) (L) we can extract 7 as follows: by its very definition, 7 is the critical 
exponent controlling the leading non- analytical behaviour of the partition function 
Z^ eu \K) as a function of the cosmological constant A. From definition (|l|) and 
the arguments given above, differentiating Z^ eu \K) twice with respect to A leads 
to the partition function where we sum over closed surfaces of spherical topology 
with two marked points. In ( |5Tf ) we have already marked one point on the disc by 
differentiating with respect to A. All that remains to be done in order to create 
closed surfaces with two marked points is to divide W^ u \l) by L to remove a 
factor proportional to the boundary length L (which was originally introduced for 
combinatorial convenience), and to contract the boundary loop L to zero. This 
gives 

d wt\L) J_ VX 
OA L L3/2 £1/2 + " ' • 

Expanding this expression in L in the limit L — > 0, and extracting the first non- 
analytic power of A leads to the critical exponent 

1 = ~\- (63) 

Following the analogous procedure for the disc amplitude ( |58| ) leads to 7 = 1/2. 

However, the physical interpretation of 7 in our model cannot be the same as in 
Liouville gravity. Firstly, the standard interpretation of 7 = 1/2 would be that we 
are dealing with objects with the fractal structure of so-called branched polymers 



20], which is obviously not the case. We also remind the reader that apart from 



characterizing the leading singularity in the partition function, in Euclidean quan- 



tum gravity the exponent 7 governs the rate of baby universe creation [21, ^2], FJ3 
This is implicit in the expansion fl62|) in that v^A is multiplied by a divergent power 
of L. This singularity reflects the proliferation of baby universes at the cut-off scale. 
By contrast, the analogous term in the expansion 

M '-vJt... (64) 

Lj Lj 
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for (p8|) is non-singular, in correspondence with the fact that by construction our 
model contains no baby universes. 

The proliferation of baby universes is closely related to the intrinsic Hausdorff 
dimension dn = 4 in Euclidean 2d quantum gravity, as decribed by the relation 
(0). As a consequence of this fractal dimensionality, the geodesic distance R, or 
equivalently the "time T" in the amplitude G£ (Li, L 2 ;T) has anomalous length 
dimension [L] 1 / 2 . This also implies that for large T the average "spatial" volume 
of a slice at some intermediate time will have an anomalous dimension. In fact, in 
Euclidean 2d quantum gravity we have 

Gt u \L u L 2 -T)(xe- nT for T — oo, (65) 

from which one can calculate the average two-dimensional volume V(T) in the en- 
semble of universes with two boundaries separated by a geodesic distance T, 

(V(T)) = 4r G [ l u \L^ L 2 ; T) oc 4tt- (66) 

For large T we therefore expect the average spatial volume L space at intermediate 
T's to behave like 

\L space) = ^ OC j^. (67) 



In the present model, according to ([44]), the amplitude behaves for large T like 

GA^L^cxe"^, (68) 

which simply means that the dimension of T in this case is [L] . We therefore obtain 
instead of ( j67|) the dependence 

(Lspace) OC -y=. (69) 

v A 



This reflects the fact that the quantum space-time of our model does not have an 
anomalous fractal dimension, and thus differs drastically from the average space- 
time in the usual two-dimensional Euclidean quantum gravity. 

It is possible to calculate explicitly (L space (T )) for the spatial volume at time 
T < T, even at the discretized level. The details are given in the appendix. 



6 Topology changes 

In our non-perturbative regularization of 2d quantum gravity we have so far not 
included the possibility of topology changes of space. At the same time we found 
disagreement with the theory of Euclidean two-dimensional quantum gravity, whose 
properties we summarized in the introduction. We will now show that if one allows 
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for spatial topology changes, one is led in an essentially unambiguous manner to the 
theory of two-dimensional quantum gravity, as defined by dynamical triangulations 
or Liouville theory. 

By a topology change of space in our Lorentzian setting we have in mind the 
following: a baby universe may branch off at some time T and develop in the future, 
where it will eventually disappear in the vacuum, but it is not allowed to rejoin 
the "parent" universe and thus change the overall topology of the two-dimensional 
manifold. This is a restriction we impose to be able to compare with the analogous 
calculation in usual 2d Euclidean quantum gravity. 

It is well-known that such a branching leads to additional complications, com- 
pared with the Euclidean situation, in the sense that, in general, no continuum 
Lorentzian metrics which are smooth and non-degenerate everywhere can be de- 



fined on such space-times (see, for example, [24] and references therein). These 
considerations do not affect the cosmological term in the action, but lead poten- 
tially to contributions from the Einstein-Hilbert term at the singular points where 
a branching or pinching occurs. 

We have so far ignored the curvature term in the action since it gives merely 
a constant contribution in the absence of topology change. We will continue to do 



so in the slightly generalized setting just introduced. The continuum results of f24 
suggest that the contributions from the two singular points associated with each 
branching of a baby universe (one at the branching point and one at the tip of 
the baby universe where it contracts to a point) cancel in the action. The physical 
geometry of these configurations may seem slightly contrived, but they may well be 
important in the quantum theory of gravity and deserve further study. However, for 
the moment our main motivation for introducing them is to make contact with the 
usual non-perturbative Euclidean path-integral results. 

We will use the rest of this section to demonstrate the following: once we allow 
for spatial topology changes, 

(1) this process completely dominates and changes the critical behaviour of the 
discretized theory, and 

(2) the disc amplitude W\(L) (the Hartle-Hawking wave function) is uniquely 
determined, almost without any calculations. 

Our starting point will be the discretized model introduced in sec. |[ Its disc ampli- 
tudes will be denoted by w^(l,g) and w^ b '(x,g), where the superscript ^ indicates 
the "bare" model without spatial topology changes. Similarly, the transfer matri- 
ces will be labelled by G^(li,l2',t = 1) and Gx(h,h',t = 1), and the continuum 
amplitudes by w[ b) (L) } wf\x). 

There are a number of ways to implement the creation of baby universes, some 
more natural than others, but they all agree in the continuum limit, as will be 
clear from the general arguments provided below. We mention just two ways of 
implementing such a change. The first is a simple generalization of the forward step 
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Figure 3: A "baby universe" branches off locally in one time-step. 




Figure 4: An alternative representation of the process in fig. |3|: a "baby universe" 
is created by a global pinching. 

we have used in the original model, where each vertex at time t could connect to n 
vertices at time t + 1. We now allow in addition that these sets of n vertices (for 
n > 2) may form a baby universe with closed spatial topology S 1 , branching off 
from the rest. The process is illustrated in fig. An alternative and technically 
somewhat simpler way to implement the topology change is shown in fig. stepping 
forward from t to t + 1 from a loop of length l± we create a baby universe of length 
I < lx by pinching it off non-locally from the main branch. We have checked that 
the continuum limit is the same in both cases. For simplicity we only present the 
derivation in the latter case. 

Accounting for the new possibilities of evolution in each step according to fig. |], 
the new and old transfer matrices are related by 

h-i 

G x (h, l 2 ; 1) = Gf (/!, l 2 - 1) + IMk-l, g) G®{1, l 2 ; 1). (70) 

i=i 

The factor li in the sum comes from the fact that the "pinching" shown in fig. |] 
can take place at any of the l\ vertices. As before, the new transfer matrix leads to 
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new amplitudes G x (li,l 2 ;t), satisfying 



and in particular 



G\(h, W, ti + t 2 ) — X! G x (h, I] t%)G\(l, l 2 ] t 2 ), 
1 



G x (h,l 2 ;t) = Y,G x (h,l;l) G x (l,l 2 ;t-1) 
1 



(71) 



(72) 



Performing a (discrete) Laplace transformation of eq. ([F^) leads to 



G(x,y;g;t) 



dz 
2iri z 



G[ b \x, z- 1 ; 1) + x^- (w(x; g)G x b \x, z~\ 1)) 



dx 



or, using the explicit form of the transfer matrix G b \x, z; 1), formula ( |17|) 



G(z,y;g;t-1), 
(73) 



G{x,y;g;t) 



dw{x,g) , J 

1 + z q h xw(x,g) — 

ox ox 



gx 



G 



1 — gx ^l—gx 



9 



(74) 



At this point neither the disc amplitude w(x,g) nor the amplitude G(x,y; g;t) are 
known. We will now show that they are uniquely determined if we assume that 
the boundary length scales canonically with the lattice spacing, L = a I, imply- 
ing a renormalized boundary cosmological constant X with the dimension of mass, 
x = x c (l — aX). In addition we assume that the dimension of the renormalized 
cosmological constant A is canonical too, g = g c (l — |Aa 2 ). Somewhat related 
arguments have been presented in different settings in 
It follows from relation ([711) that we need 



G x (h,l 2 ,t)^aG A (L u L 2 ;T). 



(75) 



It is important for the following discussion that G x (li,l 2 ;t) cannot contain a non- 
scaling part since from first principles (subadditivity) it has to decay exponentially 
in t. By a Laplace transformation, using x = x c (l — aX) in the scaling limit, we 
thus conclude that 

G x (x,l 2 ,t)^G A (X,L 2 ,T), (76) 
and further, by a Laplace transformation in L 2 , 

a->0 



G x (x,y;t) a -^ia- 1 G A (X,Y;T). 



(77) 



We will now show that the scaling of w(x, g) is quite restricted too. The starting 
point is a combinatorial identity which the disc amplitude has to satisfy. The argu- 
ments are valid both for the disc amplitude in Euclidean quantum gravity and the 
disc amplitude we have introduced for our model in (58). The discretized version of 
formula (l58f) is 



w^ b \x, g) := £ GW(x, l 2 = 1; g- 1) = G® (x, l 2 = l;g). 



(7E 
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It follows from eq. ([76]) that 



w {b \x,g) ->tr 1 Wl 6) (X). (79) 

This scaling is indeed very different from the scaling of the disc amplitude in Eu- 
clidean 2d gravity where one has 

w^(x,g) = wff>(x,g) + a^W^{X). (80) 

In relation (jSOj) , w^\x,g) is the non-scaling, analytic part of w^ eu \x,g), and 
is given by (|61|). We will assume the general form 

w(x, g) = w ns (x, g) + a v W A (X) + less singular terms (81) 

for the disc amplitude. In the case rj < the first term is considered absent (or 
irrelevant). However, if rj > a term like w ns will generically be present, since any 
slight redefinition of coupling constants of the model will produce such a term if it 
was not there from the beginning. 

We will introduce an explicit mark in the bulk of w(x,g) by differentiating with 
respect to g. This leads to the combinatorial identity 

9 ^P^- = E E G(x, I; g; t) I w(l, g), (82) 
d 9 ti 

or, after the usual Laplace transform, 

dw{x,g) r dz _ t dw(z,g) 

9 ^g— = ^f2^- Z G{X > Z ]9]t) ^z— (83) 

The situation is illustrated in fig. ^. A given mark has a distance t (T in the 
continuum) to the entrance loop. In the figure we have drawn all points which have 
the same distance to the entrance loop and which form a connected loop. In the bare 
model these are all the points at distance t. In the case where baby universes are 
allowed (which we have not included in the figure), there can be many disconnected 
loops at the same distance. Let us assume a general scaling 

T = a £ t, e>0, (84) 

for the time variable T in the continuum limit. Above we saw that the bare model 
without baby universe creation corresponded to e — 1. With the generalization (|84j) 
we account for the fact that by allowing for baby universes we have introduced an 
explicit asymmetry between the time- and space-directions. 
Inserting (|8T|) and ( jS4l) into eq. ( |8~3D we obtain 

dw ns n n _ 2 dW A (X) 1 



dg 



dw ns 1 dW A {Z) 



-aP 



dz z c dZ 



(85) 

where (x,g) = (x c ,g c ) in the non-singular part. 

From eq. (|85|) and the requirement e > it follows that the only consistent 
choices for r] are 
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T 




Figure 5: Marking a vertex in the bulk of W A (X). The mark has a distance T from 
the boundary loop, which itself has one marked vertex. 



1. rj < 0, i.e. 



1 dW A (Z) 



Z r 8Z 



(86) 



in which case we get e = 1; and 
2. 1 < r] < 2. Here formula (|85|) splits into the two equations 



and 



<9uv, 



dW A pO 1 dw % 



dA 



2a £ dz 



JdT J dZG A (X,-Z;T), (87) 



dg 



7^-1 



9=9c 



JdT J dZ G A (X, -Z; T) 



1 dW A {Z) 
7 r ~dZ~' 



We are led to the conclusion that e — 1/2 and = 3/2, which are precisely 
the values found in Euclidean 2d gravity. Let us further remark that eq. 
in this case becomes 



dW A (X) 
OA 



const. G A (X, L 2 = 0), 



which differs from (|58|), but agrees with d59|) . Finally, eq. (|88D becomes 

dW A (Z) 



JdT JdZ G A (X,-Z;T) 



dZ 



const. 



?9) 



(90) 



which will be satisfied automatically if 77 = 3/2 and e = 1/2, as we will show 
below. 
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At first sight it appears surprising that we cannot have values 77 > 2, since it is 
known that the so-called multi-critical matrix models |27] have r] = m — 1/2, m > 2 
(m = 2 corresponds to pure gravity). However, in these situations a generic coupling 
constant g does not correspond to a cosmological constant, and differentiation with 
respect to g has a different meaning. The scaling in these theories is therefore 
different (they describe non-unitary matter coupled to 2d quantum gravity and 
have negative-dimensional operators which dominate over the cosmological constant 
term). 

We will now analyze a possible scaling limit of (|73|), assuming the canonical 
scaling x = x c (l — aX) and g — g c (l — |Aa 2 ). In order that the equation have a 
scaling limit at all, x c , g c and w ns (x c ,g c ) must satisfy two relations which can be 
determined straightforwardly from (|74|). The remaining continuum equation reads 

a £ ^G A (X,F;T) = -a A [(X 2 - A)G A (X, Y; T) 



-dP- 



-i_d_ 
dX 



W A (X)G A (X,Y;T) 



(91) 



The first term on the right-hand side of eq. (0) is precisely the one we have already 
encountered in our original model, while the second term is due to the creation 
of baby universes. Clearly the case rj < (in fact i] < 1) is inconsistent with 
the presence of the second term, i.e. the creation of baby universes. However, since 
7] < 2, the last term on the right-hand side of ([)!]) will always dominate over the first 
term. Once we allow for the creation of baby universes, this process will completely 
dominate the continuum limit. In addition we get e = r\ — 1, in agreement with Q5B|). 
It follows that 7] > 1 and we conclude that e = 1/2, 77 = 3/2 are the only possible 
scaling exponents if we allow for the creation of baby universes. These are precisely 
the scaling exponents obtained from two-dimensional Euclidean gravity in terms of 
dynamical triangulations, as we have already remarked. The topology changes of 
space have induced an anomalous dimension for T. If the second term on the right- 
hand side of ( |9~ID had been absent, this would have led to e — 1, and the time T 
scaling in the same way as the spatial length L. 

In summary, in the case (r],e) = (3/2,1/2) eq. (|9l|) leads to the continuum 
equation 

d „ _ d 



— G A (X,Y;T) 



OX 



W A (X)G A (X,Y;T) 



(92) 



which, combined with eq. (0), determines the continuum disc amplitude W A (X). 
Integrating (^) with respect to T and using that G A (Li, L 2 ; T = 0) = 6(Li — L 2 ), 
i.e. 

G A (X,L 2 = 0;T = 0) = 1, (93) 



we obtain 



- 1 



_d_ 

dX 



W A (X)-^W A (X) 



(94) 
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Since W A (X) has length dimension -3/2, i.e. Wl(X) = X 3 F(y/X/X), the general 
solution must be of the form 



W A (X) = V-2AX + b 2 X 3 + c 2 A 3 / 2 . (95) 

From the very origin of W\(X) as the Laplace transform of a disc amplitude W\(L) 
which is bounded, it follows that W\(X) has no singularities or cuts for ReX > 0. 
This requirement fixes the constants b, c in (|95|) such that 



s/2 



(96) 



where the constant b is determined by the model-dependent constant in (|89|) . This 
expression for the disc amplitude agrees after a rescaling of the cosmological con- 
stant with W A (X) from 2d Euclidean quantum gravity. With W\{X) substituted 
into fl9"2"|), the resulting equation is familiar from the usual theory of 2d Euclidean 
quantum gravity, where it has been derived in various ways [p8|, |25| , [261 , with T 
playing the role of geodesic distance between the initial and final loop. In particular, 
the intrinsic Hausdorff dimension is dn = 4 as soon as we allow for baby universe 
creation. 

Let us finally comment on the difference between the equations for the amplitudes 
(|56])-(P%|) for (r],e) = (—1,1) and (f],e) = (3/2,1/2) respectively. In the first case 
there are no baby universes and eq. fl86Q entails that only macroscopic loops at a 
distance T from the entrance loop are important (as illustrated by fig. |5[). On the 
other hand, the term dW\(Z)/dZ which describes the presence of these macroscopic 
loops is absent in eq. (^). This is consistent with eq. fl89[) , which shows explicitly 
that the length of the upper loop in fig. |^ remains at the cut-off scale, i.e. it never 
becomes macroscopic. This agrees with the dominance of baby universes: at any 
point in space-time the probability for creating a little "tip" of the size of the cut-off 
scale will dominate. At the same time the right-hand side of eq. (|86|), i.e. eq. (|90D 
will play no role in the case 1 < 77 < 2, being simply equal to a constant. This latter 
property is satisfied automatically, as can be seen by using an equation analogous to 
([92]) for the exit instead of the entrance loop. Thus eq. (|90|) becomes proportional 
to 

/°°dT^ G(X,L 2 = 0;T) = const., (97) 
Jo 01 

proving our previous assertion. 



7 Conclusions 

We have tried to construct a non-perturbative model of 2d Lorentzian quantum 
gravity for universes with cylindrical topology and compact space-like slices. The 
regularization of the model was performed in the spirit of dynamical triangulations, 
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with fixed edge lengths for each triangle. To encode the light-cone structure of the 
Lorentzian geometries, we have restricted the sum over states to metric configu- 
rations with a discrete causal structure. A class of such causal triangulations was 
constructed by varying the connectivity according to the rules laid down in sec. |[ 
If we regard the edge lengths as invariant geodesic distances on the piecewise 
linear manifold corresponding to the triangulation, each triangulation defines a ge- 
ometry, i.e. an equivalence class of metrics, constructed according to Regge's pre- 



scription [29]. In this way the class of triangulations we are considering constitutes 



a grid in the set of all geometries allowing for a causal structure. As in the case of 
dynamical triangulations, we conjecture that this subset becomes uniformly dense 
at a critical point, where a continuum limit can be taken. As a consequence, neither 
a gauge-fixing nor a Faddeev-Popov determinant are needed. We should stress that 
this is an assumption, as it is in the framework of dynamical triangulations. In 
this latter approach it is corroborated by the fact that its results agree with the 
continuum theory, whenever they can be compared. Also in the present model the 
assumption seems justified, since we obtain agreement with continuum calculations, 
and, when allowing the creation of baby universes, with the dynamical triangulation 
results. 

A continuum limit of our model exists if we permit an analytic continuation in the 
coupling constants. Our results then agree with the formal continuum calculations 
performed in the so-called proper-time gauge |18[]. However, they disagree with 
those of dynamical triangulations, even if we stay "Euclidean" , i.e. do not continue 
the cosmological constant A — > —iA back to its original "Lorentzian" value. For 
example, we obtain du = 2 for the intrinsic Hausdorff dimension, which shows 
that the typical geometry of a configuration entering in the path integral is much 
smoother than in the usual Euclidean version of dynamical triangulations. We have 
therefore shown that there exists a consistent and non-trivial theory of pure two- 
dimensional quantum gravity, which can be defined as the continuum limit of a 
discrete path integral, and which does not lie in the universality class of the usual 
Liouville gravity. However, as discussed in sec. |], once we admit baby universe 
creation such that the spatial topology can change, we are forced back into the 
universality class of field theories represented by Euclidean 2d quantum gravity. 

The comparison between the Lorentzian and Euclidean sectors of our model is 
not without subtleties. One can define quantities in the Lorentzian sector that rely 
on the distinction between space- and time-like directions, and which do not possess 
an immediate analogue in the Euclidean theory. A related observation is that from 
a purely Lorentzian viewpoint, our extension to configurations with branching baby 
universes is not forced upon us from first principles. Apart from introducing points 
where a branch of the universe may disappear into nothing, the two-loop propagator 
G\(Li, L 2 ;T) has a peculiar "acausal" property. Although the causal structures on 
the individual histories are still well defined in the case of topology change, the baby 
universes that contribute to the amplitude G\(Li, L 2 ; T) consist entirely of vertices 
that do not lie in the past of any part of the central universe, (except for a single 
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point of zero volume). Their branches can even extend to times T' > T. 

We could of course equally well have worked in a "time-reflected" picture where 
the baby universes are branches coming from the past which join the central universe 
at some later stage. This would have avoided the problem with the propagator, but 
we would have had to allow for the spontaneous creation of baby universes. We 
do not think that the distinction matters for our present purposes, or that these 
features are a reason for serious concern at this stage. Our main motivation for the 
extension of our model was to understand which part of the construction needs to 
be modified in order to make contact with the usual Euclidean results. Since the 
Lorentzian configurations possess an additional, causal structure, one has to expect 
subtleties of the kind just mentioned. 

Taking the possibility of topology changes more seriously would seem to require 
a third quantized version of gravity (equivalently, a string field theory), in order to 
deal in a consistent way with the creation and annihilation of universes of length 
L. In the context of (Euclidean) non-critical string field theory, much progress has 
been made, in particular for c = 0, that is, pure 2d Euclidean quantum gravity 



2q , |26| , pOf . However, its relation with a Lorentzian theory is totally unclear. Is 
any such formulation consistent with a concept of causality? Can one implement a 
restriction to causal structures in the Lorentzian path integral? Will the pathologies 
of Lorentzian metrics at points of topology change play an important role in the 
path integral? (Unlike in the special case of baby universes, one does not in general 
expect a cancellation among the action contributions coming from such curvature 
singularities.) Maybe one will be forced to include geometries of both Euclidean 
and Lorentzian signature in the path integral, as was suggested in |2|. This work 
was motivated by some subtleties associated with the analytic continuation in the 
proper time, similar to those we have encountered in our transfer-matrix approach. 

One may summarize our results as follows: the problem of defining a theory of 
two-dimensional quantum gravity as the continuum limit of a discrete path integral 
does not have a unique answer. Addressing the problem in a Euclidean setting, 
one obtains Liouville gravity. In a Lorentzian framework, summing only over met- 
rics with a causal structure, the simplest consistent model leads to a different, in- 
equivalent theory (even modulo any analytic continuation). If one wants to obtain 
agreement between the two, one must extend the Lorentzian model by allowing for 
spatial topology changes. However, from a purely Lorentzian point of view it is 
unclear whether or not such configurations should be included in the path integral. 

What can we conclude from the present work for the case of real four-dimensional 
Lorentzian gravity? One clearly can formulate similar causal restrictions on dynam- 
ically triangulated four-geometries, but we do not know whether this prescription 
alone would make the path integral sufficiently Lorentzian. The presence of the 
Einstein-Hilbert term makes the issue of analytic continuation more complicated. 
However, if our results in two dimensions are anything to go by, summing over ge- 
ometries with a causal structure could in principle lead to a drastic change of the 
results. It may be worthwhile to test this idea by numerical simulations, which are 
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well-developped for the case of Euclidean geometries (see, for instance, |3T| and |32 
for most recent developments). 
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Appendix 

In this appendix we calculate the average size of the spatial universe at time Ti, 
when the total time is given by T\ + T 2 . We will show that in the continuum 
limit the spatial slices are genuinely extended with a spatial volume ~ 1/ y/X. One 
can use various boundary conditions. The computation is simplest when one fixes 
the boundary cosmological constants X, Y and lets the boundary lengths fluctuate 
according to the distribution dictated by the choice of the boundary cosmological 
constants. The calculation can be done even at a discretized level. It is helpful to 
use the properties 

A tl A t2 - B h B t2 = (1 - F 2 )A tl+t2 , B tl B t , 2 - C tl C t2 = (1 - F 2 )C tl+ t 2 (98) 



Bt x A t2 — C tl B t2 = A tl B t , 2 — B tl C t2 — (1 — F )B tl+t2 (99) 

of the coefficients A t , B t and Ct, which follow directly from the fact that G(x, y; g; t), 
parametrized as in ([22] ) must satisfy the fundamental composition law ([13]). 

For purely aesthetic reasons, we will consider the symmetric situation where both 
the entrance and exit loops are marked. We calculate the average size (l(ti)) xy of a 
spatial universe at time-slice t\ as 

QfaKv = 7X277 7— TT f T— G{x,z-\g-M) z—G^(z,y;g;t 2 ), (100) 

,y GW {x, y; g; ti + i 2 ) J 2m z dz 

where the contour encloses the poles in the z-plane of G(x, z~ x \ g; ti) but not of 
G^(z, y; g; t?). From (|27|) and (|98|)-(|99|) the result of the contour integral is 

(](f)) 1 (A tl - B tl x)(A t2 - B t2 y) + (B h - C tl x)(B t2 - C t2 y) 

{{l))x ' y l-F* A h+t2 - B h+t2 (x + y) + C h+t2 xy ' { ] 

If t\ and t 2 go to infinity and a to zero in such a way that T% and T 2 stay finite, at 
the critical point g c = \, but with x and y not going to 1 (the case of microscopic 
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boundaries), we have 



1 (1 _ e -2VAT l)(1 _ e -2VAT 2 



aVA 



1 _ e -2 v / A(Ti+T 2 ) 



(102) 



This shows that (L(Ti)) = (al(ti)) xy possesses a continuum limit, independent of 
the size of the microscopic boundaries and independent of T\ and T 2 if they are 
sufficiently large. 

From ( |101| ) we can directly find the expression {L{T/)) XY in the case of macro- 
scopic loops: 



(L{Ti)) 



[\/X + X) + e" 2v/XTl (VA - X)] Uy/X + Y) + e - 2 ^ (y/X - Y) 



xy 



y/X (X + \/X)(Y + y/X) - e - 2 VA(T 1+ T 2 )(^X -X)(y/X- Y) 



(103) 

The same expression could have been obtained starting directly from the continuum 
expression (|I0"D or equivalently from 



G ( 1\X,Y;T) = 
and by calculating 

(m))x. Y = 



4A e~ 2 ^ T 



{X + y/X) (Y + y/X) — e- 2 ^ T (y/X - X)(y/X - Y) 



2 • 



(104) 



1 



Gf{X,Y-T x + T 2 ) 
1 

GfiX^Y-^ + Tz) 



IOO 

-ioo 
too 



d 



dZ G A (X, -Z; T x )-^Gf\z, Y; T 2 ) 



f °° ciZ Gf(X, -Z; T^cf } (Z, F; T 2 )(105) 

J —ioo 



Let us finally compute (L(Ti)) L La with the boundary lengths L\ and L 2 kept 
fixed. In this case there is no problem with the marking of boundary loops since the 
factors of L x and L 2 cancel in the normalization. Thus 



1 



roc 

/ dL Ga(Li, L; T\) L G\(L, L 2 ; T 2 ) (106) 
Jo 



M G A (L U L 2 ]T 1 +T 2 ) Jo 
Rather surprisingly, the integration can be performed explicitly using 



dxx 3 e— 2 / 1 (^)/ 1 ( 7 x) = ± 1 e^[{^ + 1 2 )h(^L) +2/3 7 / (^)], (107) 



2a- 



and we obtain (before normalization), with the notation Si = sinh y/XTi and C{ 
cosh a/ATj, i = 1, 2, 



/ r /'T' \ \unnorm 



,i 2 



V^gg VA(§iL 1+ f L 2 ) 

r o3 e e 
-^2>Jl+ 2 



X 



(L 1 ^ 2 2 + J L 2 5' 1 2 )J 1 ( 



2/AL 1 L 2 



S 



1+2 



+ 2SiS 2 \JT^L 2 Iq(^ 



2/M x L r 

Sl+2 



(10* 
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Dividing by the normalization factor we obtain 

1 



D l+2 
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